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Abstract 

By using Malliavin calculus and multiple Wiener-Ito integrals, we study the existence 
and the regularity of stochastic currents defined as Skorohod (divergence) integrals with 
respect to the Brownian motion and to the fractional Brownian motion. We consider 
also the multidimensional multiparameter case and we compare the regularity of the 
current as a distribution in negative Sobolev spaces with its regularity in Watanabe 
space. 
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1 Introduction 

The concept of current is proper of geometric measure theory. The simplest example is the 
functional 

Jo 

defined over the set of all smooth compact support vector fields ip : M. d — > R rf , with 7 : 
[0, T] — > M. d being a rectifiable curve. This functional defines a vector valued distribution. 
Let us denote it by 

5{x- 1 {t))i(t)dt. (1) 



I 

Jo 
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We address for instance to the books [3J, [23J, [TU], [TJ, for definitions, results and applica- 
tions. 

The stochastic analog of 1-currents is a natural concept, where the deterministic 
curve (l (t))te[oT] * s re pl ace d by a stochastic process (X t ) te ^ 0T ^ in R. d and the stochastic 
integral must be properly interpreted. Several works deal with random currents, see for 
instance [EJ, M, [22], [EJ, [E], [IT], [EJ 0], [EJ, [6]. 

The difference between classical integration theory and random currents is the at- 
tempt to understand the latter as random distribution in the strong sense: random vari- 
ables taking values in the space of distributions. The question is not simply how to define 
Jq 1 ((p (X t ) , dX t ) R d for every given test functions, but when this operation defines, for almost 
every realization of the process X, a continuous functional on some space of test functions. 
This is related to (a particular aspect of) T. Lyons theory of rough paths: a random current 
of the previous form is a concept of pathwise stochastic integration. The degree of regularity 
of the random analog of expression ([!]) is a fundamental issue. 

In this work we study the mapping defined as 

f(s) = / 5(x-B s )dB s , x£R d ,T>0. (2) 

J[0,T] N 

Here the integrator process is a d- dimensional Wiener process with multidimensional time 
parameter or a ci-dimensional fractional Brownian motion. 

A first direction of investigation is the regularity with respect to the variable x 6 M. d 
of the mapping given by ([2]) in the (deterministic) Sobolev spaces H~ r (W 1 ;W i ). This is the 
dual space to H r (M. d , M. d ) , or equivalently the space of all vector valued distributions <p such 

that (! + M ») - W W* to < oo, ? be ing t he Fourier tr a»s form of v . So,ne JL on 

this have been recently obtained by [4] , [5] or [6] using different techniques of the stochastic 
calculus. We propose here a new approach to this problem. The main tool is constituted by 
the Malliavin calculus based on the Wiener-Ito chaos expansion. Actually the delta Dirac 
function 5(x — B s ) can be understood as a distribution in the Watanabe sense (see [23] ) 
and it has been the object of study by several authors, as [21], pT] or [12] in the Brownian 
motion case, or [I], [2] or [8] in the fractional Brownian motion case. In particular, it is 
possible to obtain the decomposition of the delta Dirac function into an orthogonal sum of 
multiple Wiener-Ito integrals and as a consequence it is easy to get the chaos expansion of 
the integral ([2]) since the divergence integral acts as a creation operator on Fock spaces (see 
Section 2). As it will be seen, the chaos expansion obtained will be useful to obtain the 
regularity properties of the mapping £. Another advantage of this method is the fact that 
it can be relatively easily extended to multidimensional settings. 

Besides the estimation of the Sobolev regularity with respect to x, we are also 
interested to study the regularity with respect to to, that is, as a functional in the Watanabe 
sense, of the integral ([2]). Our interest comes from the following observation. Consider 
N = d = 1. It is well-known (see Nualart and Vives [21] ) that for fixed i£lwe have 

5(x -B s )e IDT"' 2 for any a > - 
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where B _a ' 2 are the Watanabe (or Sobolev-Watanabe) spaces introduced in Section 2. On 
the other side, more or less surprisedly the same order of regularity holds with respect to 
x; for fixed u, the mapping g[x) = 8{x — B s (u)) belongs to the negative Sobolev space 
H- r (R;R) for every r > \. Indeed, since the Fourier transform of g is g(x) = e lxBa we 
have 

Mh-<-(r ; r) = / \g(x)\ 2 (l + x 2 )~ r dx 

and this is finite if and only if r > \ . 

One can ask the question if this similarity of the order of regularity in the determin- 
istic and stochastic Sobolev space still holds for the functional £ defined by[2j and actually 
in the case of dimension d = 1 the above property still holds: we have the same regularity 
of £ both with respect to x and with respect to uj. 

One can moreover ask if it holds for other Gaussian processes. The answer to this 
question is negative, since we show that (actually, this has also been proved in [6] but 
using another integral with respect to fBm) in the fractional Brownian motion case the 
mapping ([2]) belongs to ff~ r (IR; M) for any r > ^ — ^ (as a function in x) and to B~ a ' 2 
with a > | — (as a function on uj). 

We organized our paper as follows. Section 2 contains some preliminaries on Malli- 
avin calculus and multiple Wiener-Ito integrals. Section 3 contains a discussion about 
random distribution where we unify the definition of the quantity 5(x — X s (uj)) (X is a 
Gaussian process on M, d ) which in principle can be understood as a distribution with re- 
spect to x and also as a distribution in the Watanabe sense when it is regarded as o function 
of uj. In Section 4 we study the existence and the regularity of the stochastic currents driven 
by a N parameter Brownian motion in M. d while in Section 5 concerning the same prob- 
lem when the driving process is the fractional Brownian motion. In Section 6 we give the 
regularity of the integral ([2]) with respect to uj and we compare it with its regularity in x. 

2 Preliminaries 

Here we describe the elements from stochastic analysis that we will need in the paper. 
Consider TL a real separable Hilbert space and (B(p), (p £ TL) an isonormal Gaussian process 
on a probability space (f2, A, P), that is a centered Gaussian family of random variables such 
that E (B((p)B(ip)) = (ip,ip)i-[. Denote by I n the multiple stochastic integral with respect 
to B (see |20|). This I n is actually an isometry between the Hilbert space 7i &n (symmetric 
tensor product) equipped with the scaled norm -p=|| • and the Wiener chaos of order 

n which is defined as the closed linear span of the random variables H n (B(ip)) where tp S 
TC, \\<p\\h = 1 and H n is the Hermite polynomial of degree n > 1 
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The isometry of multiple integrals can be written as: for m, n positive integers, 



E(I n (f)I m (g)) = n\{f,g) n ®n if m = n, 

E(I n (f)I m (g)) = ifm/n. (3) 

It also holds that 

In(f)=In{f) 

where / denotes the symmetrization of / defined by f(xi,. . . , x x ) = ^ Y,aeS n f( x a(l),- ■ ■ , %*{n)) 

We recall that any square integrable random variable which is measurable with 
respect to the cr-algebra generated by B can be expanded into an orthogonal sum of multiple 
stochastic integrals 

F=Y,Ufn) (4) 

n>0 

where f n E 7i® n are (uniquely determined) symmetric functions and Jo(/o) = E [F]. 
Let L be the Ornstein-Uhlenbeck operator 



LF = -J2 nI n(fn) 



n>0 

if F is given by 

For p > 1 and a£lwe introduce the Sobolev-Watanabe space D a,p as the closure 
of the set of polynomial random variables with respect to the norm 

H^lkp = ~ L)%\\ LP{n) 

where I represents the identity. In this way, a random variable F as in 

(HD belongs B a > 2 if 

and only if 

X)(l + <[IW»)llia ( n)=E( 1 + n ) an! ll/«llw»» <00 - 

n>0 n>0 

We denote by D the Malliavin derivative operator that acts on smooth functions of the 
form F = g(B(ipi), . . . , B((p n )) [g is a smooth function with compact support and ifi £ 7i) 

DF = Y J -^{B(vi)T--,B^ n ))v l . 

i=l 1 

The operator D is continuous from B Q ' p into B Q_1 ' P (7i) . The adjoint of D is denoted by 5 
and is called the divergence (or Skorohod) integral. It is a continuous operator from O a ' p (TC) 
into D a_1,p . For adapted integrands, the divergence integral coincides to the classical ltd 
integral. We will use the notation 



5(u) = / u s dB s . 
Jo 
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Let u be a stochastic process having the chaotic decomposition u s = Y2 n >o (/«(*> •?)) where 
/«,(•> s) G J{® n for every s. One can prove that u G Dom 5 if and only if f n G 7^®( n+1 ) for 
every n > 0, and X)^=o -^n+l(/n) converges in L 2 (f2). In this case, 

oo oo 

<5(u) = ^W/„.) and E|5(n)| 2 = J> + 1)! ||/„||^ (n+1) . 

n=0 n=0 

In the present work we will consider divergence integral with respect to a Brownian 
motion in R d as well as with respect to a fractional Brownian motion in R . 

Throughout this paper we will denote by p s (x) the Gaussian kernel of variance s > 

2 

given by p s (x) = } e~2i, x G M and for x = (x\, . . . , Xd) G R d by p d {x) = Y\ d =1 p s (xi). 



3 Random distributions 

We study now the regularity of the stochastic integral given by ([2]). Our method is based on 
the Wiener-Ito chaos decomposition. Let X be isonormal Gaussian process with variance 
R(s,t). We will use the following decomposition of the delta Dirac function (see Nualart 
and Vives [21], Imkeller et al. [11], Eddahbi et al. [2]) into orthogonal multiple Wiener-Ito 
integrals 

5(x -X S ) = J2 R(s)^p Ris) (x)H n ( * ) I n (lgjj) (5) 
n>0 \R\ S ) 2 / 

where R(s) := R(s,s), Pr( s ) is the Gaussian kernel of variance R(s), H n is the Hermite 
polynomial of degree n and I n represents the multiple Wiener-Ito integral of degree n with 
respect to the Gaussian process X as defined in the previous section. 

A key element of the entire work is the quantity 5(x — B s ). This can be understood 
as a generalized random variable in some Sobolev-Watanbe space as well as a generalized 
function with respect to the variable x. The purpose of this section is to give an unitary 
definition of the delta Dirac as a random distribution. 

Let D be the space of smooth compact support functions on M. d and let D 1 
be its dual, the space of distributions, endowed with the usual topologies. We denote by 
(S, ip) the dual pairing between S G D' and ip G D (R ) . We say that a distribution 
S G D'(R d ) is of class £f oc (M d ), p > 1, if there is / G L} oc (R d ) such that (S,<p) = 
J K d / (x) </? (x) dx for every ip G D . In this case we also say that the distribution is 
given by a function. 

Let (p,,A,P) be a probability space, with expectation denoted by E. 

Definition 1 We call random distribution (on M. d , based on (Q,A,P)) a measurable map- 
ping uj i — > S (lo) from ($7, A) to D' with the Borel a -algebra. 

Given a random distribution S, for every tp G D the real valued function 
uj i— ► (5 (to) , tp) is measurable. The converse is also true. 
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Similarly to the deterministic case, we could say that a random distribution S is 
of class L° (ft,Lj oc (R d )) if there exists a measurable function / : !] x l d -> 1, / (a;, •) G 
(^ d ) for ^- a - e - ^ G such that (S (uj) , <p) = j Rd f (u, x) <p (x) dx for every ip £ D (R d ) . 
In this case we also say that the distribution is given by a random field. However, this 
concept is restrictive if we have to deal with true distributions. 

There is an intermediate concept, made possible by the simultaneous presence of the 
two variables uj G ft and x G M. d . One could have that the random distribution is given by a 
function with respect to the x variable, but at the price that it is distribution-valued (in the 
Sobolev-Watanabe sense) in the uj variable. To this purpose, assume that a real separable 
Hilbert space H is given, an isonormal Gaussian process (W (h) ,h 6 H) on (ft,A,P) is 
given. For p > 1 and a G M, denote by B Q ' P the Sobolev-Watanabe space of generalized 
random variables on (ft, A, P), defined in terms of Wiener chaos expansion. 

Definition 2 Given p > 1 and a < 0, a random distribution S is of class Lj oc (M. d ; B a ' p ) 
i/E [\(S, ip)\ p ] < oo for all ip e D (R d ) and there exists a function f G h\ oc (M. d ; B a >P) such 
that 

{S, ip) = <p(x)f (x) dx 

for every ip G D (IR^) . 

Let us explain the definition. The integral J Rd ip (x) f (x) dx is of Bochner type (the 
integral of the D Q ' p -valued function x i— ► ip (x) f (x)). The integral j Rd <p (x) f (x) dx is thus, 
a priori, an element of B a ' p . The random variable uj <— >■ (S (uj) , ip), due to the assumption 
E [\(S, ip)\ p ] < oo, is an element of B a ' p for all a' < 0. The definition requires that they 
coincide, as elements of B a,p , a priori. A fortiori, since uj i— > (S (uj) , ip) is an element of 
B°' p , the same must be true for J Rd ip (x) f (x) dx. Thus, among the consequences of the 
definition there is the fact that J Rd ip (x) f (x) dx G B 0,p although / (x) lives only in B Q ' P . 

The following theorem treats our main example. Given any d- dimensional random 
variable X on (ft, A, P), let S be the random distribution defined as 

(S(uj),ip):=ip(X(uj)), ipGD(R d ). 

We denote this random distribution by 

5(x-X(uj)) 

and we call it the delta-Dirac at X (uj). Denote Hermite polynomials by H n (x), multiple 
Wiener integrals by I n , Gaussian kernel of variance a 2 by p a 2 (x). 

If Hi, ■ ■ ■ ,TCd are real and separable Hilbert spaces, a d- dimensional isonormal pro- 
cess is defined on the product space (ft, A, P), ft = ft\ x . . . x ftj, A = Ai®...®Ad, P = P\0 
...®Pd as a vector with independent components ((H /1 (/ii), . . . , W d (ipd)), h\ G Hi, . . . , hd G Hd) 
where for every i = 1, . . . , d, (W l (hi), hi G Hi) is an one dimensional isonormal process on 
(fti, Ai, Pi). We denote by Ej the expectation on (fti,Ai, Pi) and by E the expectation on 
(ft,A,P). 
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Theorem 1 Let Hi, i = 1, . . . ,d be real separable Hilbert spaces. Consider 
((W^/ii), . . . , W d (ifd)), hi G Hi, ■ ■ ■ , hd € H<i) a d- dimensional isonormal Gaussian process 
on (Q,A,P) as above. Then, for every hi G Hj, i = l,...,d, the random distribution 
5 (• — W(/t)) is of class Lj oc (M. d ; ID ' 2 ) /or some a < and f/ie associated element f of 
L] oc (R d ; B a > 2 ) zs 

/(*) = E ni^r^p^^fe)^^) 

m,...,n d >0i=l Vl i|y 

where x = {x\, . . . , xj) G M d and i^., z = l,...,d denotes the multiple Wiener-Ito integral 
with respect to the Wiener process 

We denote this B a ' 2 -valued function f (x) by 5 (x — W (h)). 

Proof: Let us consider first the case d = 1. In this case 



f( X ) = J2 \h\- n P\ h f (*) (^) In • 



n>0 

We have to prove that 

<p(W(h)) = f <p{x)f{x)dx 

JR 

for every test function (p. We have 

j^(x)f(x)dx = ^H" n / n (^)jf ^(Xjpi^ffn^dx 
n>0 ^ R 

On the other hand, by Stroock's formula we can write 

<p(W(h)) = ^- } I n (D^ l p(W(h)))=^- ] I n (h^)EU n \w(h))) 

TL . \ / TL . \ / 



n>0 n>0 



and 



EU n \W(h)j) = ^ n \x)p\ h? {x)dx = \h\ ^ n \x\h\) Pl {x)dx 

x ' JR JR 

= (-1) / ^ n - l \x\h\){pi{x)) ! dx 
Jr. 

= i-irihr 1 [ V (x\h\)(pi(x))^dx. 

JR 
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By using the recurrence formula for the Gaussian kernel p\ 

(Pi(x)) (n) = (-l) n n\ Pl (x)H n (x) 

we get 



<P(W{h)) = Y, [ <p(x\h\)H n (x) Pl (x)dx = f tp(x)f(x)dx. 



Concerning the case d > 2, note that 

d 



[ <p{x)f{x)dx = E fl(\hi\- ni ^A h f ni ) I dx^ X )p lhil ,( Xl )H ni (^ 
jRd m,...,n d >0i=i V jRd V|/l ' 

E II ( l^l" n<+1 4 {hf ni ) [ <p(\h 1 \x 1 ,...,\h d \x d )p 1 (x i )H ni (x i )dx 

n u ...,n d >0i=l V Jr " 

Let us apply Stroock 's formula in several steps. First we apply Stroock formula with 
respect to the component W l 

<p{W{h)) = ^(WHhi), W d {h d )) = £ -L4 (Ex [d^MW 1 ^), W d (h d ))] ) 



ni! 

ni>0 



ni>0 



and then with respect to W 



2 



y^iwHh!), w d {h d )) = e - 1 4(/ i r 2 )E 2 ^j-^ (w 1 ^), . . . , w^m). 

1 n 2 >0 2 ' V 1 2 / 



We will obtain 



v{W{h)) = E II^4(^ n ')E Q n[ J {W\h l ),... 1 W\h d ))) 

m,...,n d >0i=l l ' \ 1 " d / 

m,..,%>Oi=l 1 « 



and now it suffices to follow the case d = 1 by integrating by parts and using the recurrence 
formula for the Gaussian kernel. ■ 
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4 Regularity of Brownian currents with respect to x 

Let (Bt)t£[o t T\ N be a Gaussian isonormal process. For every n > 1, x G R and s = 
(si, . . . , s^) E [0, T]^ we will use the notation 

= fl(*)-S PiiW (s)ff n ( -^r ] (6) 

with R(s) = ES^ = si . . . Sjv- We will start by the following very useful calculation which 
plays an important role in the sequel. 

Lemma 1 Let a n be given by (Ejj and denote by a n {s) the Fourier transform of the function 
x — > a n (s). Then it holds that 

a x n (s) = e 2 2- (7) 

Proof: Using formula ([5]) one can prove that that Fourier transform of g(x) = 5(x — B s ) 
admits the chaos expansion 

^) = E^W j »( 1 m) ( 8 ) 

n>0 

where a^(s) denotes the Fourier transform of the function x —* a n (s). 

But on the other hand, we know that 5{x — B s ) = e~ lxBa and using the Stroock's 
formula to decompose in chaos a square integrable random variable F infinitely differentiable 
in the Malliavin sense 

F = V —I n \e(D^F) 

n>0 

where denotes the nth Malliavin derivative and by the trivial relation 

D u e~ ixB ° =-ixe- ixB °l [QA (u) 



which implies 

Dil., Un e- ixB ° = (-ixTe-^lf^m, ...,u n ) 

one obtains 



n>0 ' n>0 

Now, by putting together ([8]) and © we get j7|). ■ 
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4.1 The one-dimensional Brownian currents 

Let B be in this part a the Wiener process. Then by Q 



6(x -B S ) = J2 s-hs(x)H n (4) I n (lgjj) = ^ <( S )I„ (l 

n>0 \S2/ n>0 



[o,*]y ' 



(10) 



here I n representing the multiple Wiener-Ito integral of degree n with respect to the Wiener 
process. In this case we have 



<(*) = s 2 Ps(x)H n — . 



X 



S2 

We prove first the following result, which is already known (see [4], [6]) but the 
method is different and it can also used to the fractional Brownian motion case and to the 
multidimensional context. 

Proposition 1 Let £ be given by (Q|). For every it holds that 



E 



fa) 



T. 



As a consequence £ G H r (M, R) i/ and on/?/ if r > i. 

Proof: We can write, by using the expression of the Skorohod integral via chaos expansion 

C(x) = ^I n+1 ((<(,)!£,(•)) ~) (11) 



n>0 



where ^o£(s)lj^™](«) 1 denotes the symmetrization in n+1 variables of the function (s,ti, . . . , t r , 
a^(s) l^"j(ii, . .. ,t n ). It can be shown that 



|(x) = ^J n+1 (((at( S )l^ s] (- 



n>0 



By using ((TJ), 

fro = E 7 < 



n+l 



(— i) n x r 



n 



(— i) n x 



n>0 v v ' ' n>0 

and therefore, by the isometry of multiple integrals ([3]), 

2 

n>0 



E 



2 T 2n / 2 \ ' 

&)| =E^j)5( n + 1 ) , ll( e " Vl M(-)) 



l.r, 2 ([0,r]"+ 1 )- 



(12) 
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Since 



l&O) (h,...,t n+1 ) 



-I n +l 2 
1 \ - 

e 2 



n + 1 ^ [0,*i 



-i (gin / , r , 

fo t-] \ i' ■ ■ ■ ' *' ■ ■ ■ » 



n+1, 



where means that the variable U is missing, we obtain 

2 



E 



e(z) 



e 2 "''^i]^ 1 ' •••!*«)••• )*n+l)lli 2 ([0,T]' l + 1 ) 



V— i II V 

^ (n + 1)! ^ 

n>0 v ; i=0 

ETTA ^| jT]n ^---^(e-^l-(i 1 ,... ) t„) 



n>0 ,yu 
„2n f-T 



£ 

n>0 



.1' 



i.l 



Finally, 



dse x s s n ds 



ra>0 



-ds = r. 



■/?■! 



E|eiH-( R;R )= / (l+^)" r E £(*) 

JM. 

and this is finite if and only if r > i . 



(1 + x 2 )" r Tdx 



Remark 1 The above result gives actually a rigorous meaning of the formal calculation 



E (J* 8{x - B s )dB s ) =EJ* \S(x - B s )\ 2 ds = T. 

4.2 The multidimensional multiparameter Brownian currents 

Let us consider now the multidimensional situation. In this part we will actually treat the 
regularity of the function 

5(x - B s )dB s (13) 



/ 



where x £ R d and B = (B 1 , . . . , B d ) is a (f-dimensional Wiener sheet with parameter 
t G [0, T]^. Here the term (|13p represents a vector given by 



8{x - B s )dB s 

[0,T] N 

(ei(x),...,^(x)) 



If 5(x- B s )dBl, f 5(x- B s )dB d s 

\J[Q,T} N J[0,T] N , 



for every x E R rf . 

The method considered above based on the chaos expansion of the delta Dirac func- 
tion has advantage that, in contrast with the approaches in [B] or [5], it can be immediately 
extended to time parameter in M. N . Moreover, we are able to compute explicitly the Sobolev 
norm of the vector (|13|) and to obtain an "if and only if result. 
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Proposition 2 Let (B t ) te ^ ^N be a d -dimensional Wiener process with multidimensional 
parameter t G [0, T] . Then for every w G f2 the mapping x — * Jj Q T m 5(x — B s )dB s belongs 
to the negative Sobolev space H- r (R d ,R d ) if and only if r > tt. 

Proof: We need to estimate 



|2\-r 



;i + |x|T r E 



(i.r 



with 



2 



= l£i( x )| 2 + - ■ ■ l^d(^)| 2 where denotes the complex modulus of £i(x) G C. 



We can formally write (but it can also written in a rigorous manner by approximating 
the delta Dirac function by Gaussian kernels with variance e —* 0) 

d 

5(x - B.) = J] 8(xj - B{) =5 k (x- B s )8(x k - B k s ) 

3=1 

where we denoted 

d 

5 k (x-B s )= 11 5( Xj -Bi), k = l,...d. (14) 

Let us compute the kth component of the vector £(x). We will use the chaotic expansion 
for the delta Dirac function with multidimensional parameter s G [0, T] (see ([5]) 



8( Xj -Bi)= E^^)4,(C](-)J (15) 

n 3 >0 

where I n denotes the Wiener-Ito integral of order n with respect to the component B J and 
for s = (si, ...,sn) 

, 3 \{s 1 ...s N )~~. (16) 

V s l ■ ■ ■ S N / 

It holds, for every k = 1, . . . , d, by using formula (|10p 



U*)= I 5{x-B s )dB k = [ 5 k (x - B s )S(x k - B k )dB k 

5 h (x - B.) Y, °%M k (ljJJO)) ^ 



[0,T]1 



n^>0 



Since the components of the Brownian motion B are independent, the term 5 k (x — B s ) is 
viewed as a deterministic function when we integrate with respect to B k . We obtain, 

= / " = ^ I ns+1 Us k (x - B s )a*l(s)lf n s U 

J\0,T\ N ~Tr, VV 



n fc >0 
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Here ySk(x — 5 s )o^(s)l^"|(-)J denotes the symmetrization in + 1 variables of the 
function 

(s, t x , . . . , tn J = 8 k (x - B s K* (a) lgj (t ! , . . . , t nk ) . 

Let us denote by the Fourier transform of £fc(cc) = /[ 0T jjv S(x — B s )dB^. As 

in the previous section one can show that (that is, the Fourier transform with respect to x 
"goes inside" the stochastic integral ) 



Ux) = ^ 4 +1 - 5 s )a^( S )l^(-) 

n fe >0 



where 

is the Fourier transform of 



S k {x -B s )al k (s 



(17) 



Now clearly 
and by © 



Bx = (x 1 ,...,x d )^5 k (x-B s )a x n k k (s)= I ]J S(xj — B J S ) J a^(s). 



H)"^ E ( C -Mi\ _ HT fc 4 fc „4 



e 2 



with |s| = si . . . sat if s = (si, . . . , sn)- Thus relation (fT7|) becomes 

'*-E( e - fa * B . fc )lgJ*(.) 



n fc >0 



Taking the modulus in C of the above expression, we get 



COS 



+ 



X x * B 

K l=l,l-Lk 
d 



{-i) nk xl k _4W 



n k >0 



sin 



(18) 
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and using the isometry of multiple stochastic integrals ([3D and the independence of compo- 
nents we get as in the proof of Proposition 1, (we use the notation |s| = s\ . . . sjv) 



E 



60) 



E 



2n t 



X, 



n k >0 



COS 



2nu 



77. > 



sin 



n k >0 



2n k 



E 



E x i B s 

d 

E x ^ 



e 2 1 [0,4('-"lL 2 ([0,T] JV ("fe+ 1 )) 



e i [0,s] UllL2([ 0jT ]JV(™ fc +l)) 



This implies 



E 



[0,T]t 



E 

, n fc >0 



ds 



ds = T N . 



[0,7? 



dx 



l + |x| 2 )- r E( Ci(x) 



N 



[ + k 

7R d 



dx 



2\-r 



dx 



which is finite if and only if 2r > d. 



Remark 2 It is interesting to observe that the dimension N of the time parameter does 
not affect the regularity of currents. This is somehow unexpected because it is known that 
this dimension N influences the regularity of the local time of the process B which can be 
formally written as j^ T ^ N 5{x — B s )ds (see 



5 Regularity of fractional currents with respect to x 
5.1 The one-dimensional fractional Brownian currents 

We will consider in this paragraph a fractional Brownian motion (-B^)te[o,T] with Hurst 
parameter H E (i , 1). That is, B H is a centered Gaussian process starting from zero with 
covariance function 

BP{t,8) :=^{t 2H + s 2H -\t-s\ 2H ), s,t€[0,T\. 

Let us denote by TCh the canonical Hilbert space of the fractional Brownian motion which 
is defined as the closure of the linear space generated by the indicator functions {l[ 0) t],i £ 
[0, T]} with respect to the scalar product 

<l[o,t],W«ff = R H {t,s), s,te[0,T}. 
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One can construct multiple integrals with respect to B H (the underlying space 
L 2 [0,T]) is replaced by Tin) and these integrals are those who appear in the formula (fTUj) , 
We prove the following result. 

Proposition 3 Let B H be a fractional Brownian motion with H £ {\, 1) and let £ be given 
by (0). Then for every w£S), we have 



for every r > ^ — \ ■ 



Proof: Using the above computations, we will get that the Fourier transform of g{x) = 
5(x — Bf) is equal, on one hand, to 

g(x) = J2<(s 2H )I^ (ifoj(-)) 

n>0 V 

(here denotes the multiple stochastic integral with respect to the fBm B H and the 
function is defined by ([6])), and on the other hand, 



n! 

n>0 



and so 



Moreover 



We then obtain 



2 2H I A„\n 

4(s™) = e-^t^-. (19) 
n! 



n! 

n>0 



,2.2H 



2 



E ll£lltf-(IR« = E / (l + g 2)r l^)] ^ 



7 R (i + x^^ ( n + 1)! x 



n+1 ,2^2H «±1 
X (2^ e 2 1 M( U l'---:«ir",«r l +l),2 j e 2 1®"^ (^ ; . . . , Vj , . . . ; V n+1 )) Hl 
i=l j=l 

Using the fact that for regular enough functions / and g in Tin their product scalar is given 
by(seee.g. [20], Chapter 5) (/, g)n H = H(2H- 1) / T f(x)g(y)\x -y\ 2H - 2 dxdy and thus 
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for regular enough function f,g£ H® n 

(f,g) H ®n = (H(2H-l)) n / / f(u 1 ,...,u n )g(v 1 ,...,v n )T\\u i -v i \ 2H ~ 2 du 1 ..du n dv 1 ..dv T} 

H J[0,T] n J[0,T] n i=1 



we find (by duiduj we mean below du\..du n+ \dvi..dv n+ i) 

E ll?ll/f-'-(R;R) 



y R (i+ X 2)r^ („+i)! ^4^4^11^ 

f 1 Jin 

= / 7 or- Y(H(2H- l)) n+1 - rrdx 

k(l + x*y^ 1 " (n + 1)! 



n>0 

n+1 ,. ,. x 2 u 2H x' 2 v 2H ri+1 

x V / / e i~ TT |tt, - vi\ 2H ~ 2 duidv 

^J[o,r]»+ 1 i[o,T]»+i ^ 

+ Xor^E^-i)) (n + i)] 

n+1 „ „ n+1 

E / / IlK-^i 



x V / / ] I h-^ 2 "" 2 ^^ 

M ^ 3 i[0,T]»+ii[0,T]»+i ^ 

xe s 2 - 1®« ^ . . . ; ^, . . . , u n+1 )l®^ vj] (vi, , u„+i) 

:= A + B. 

Let us compute first the term A. Using the symmetry of the integrand and the fact that 

2H(2H-l)f [ \u-v\ 2H ~ 2 dudv = R(t,s), s,te[0,T] (20) 
Jo Jo 

we can write 

r 2n 

n! 



n+1 — dx 



n+1 „ ,. n+1 

xv / nwi-v^^duidv, 

2 2ff 2 2iT 

xe-^e-^l®"^, . . . , ,« n+ i)l®^ i] (w 2 , . . . Vn+ i) 
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and by integrating with respect to du2 ■ ■ ■ du n+ \. 

r 1 » x 2n 
A - L & (T+W § 2H(2H ~ v >~ri dx 

e 2 e 2 ^-^1^-" z R( u ,v) n dudv 



/• 1 f 1 f 1 1 1H 1 1H I 



2 / ^ X 2n i?(u,w) 



n! 

n>0 



//(2//-1) /' r/.r ' ., /'' [' , ^\u- ^"-l^^^dudr 



= H(2H-\) / ^- / / e 2 \u-v\ 2H - 2 dudv. 



Now, by classical Fubini, 

f T f' T f 1 i 2 |u-i,| 2H 

A = H(2H-1) / / |u - v| 2H - 2 ducfo / ^-e 2 dx 

and by using the change of variable x\u — v\ H = y in the integral with respect to , 

A = H(2H-1) f T ( T \ u -v\ 2H - 2+2Hr - H I e~ y 4(\u-v\ 2H + y 2 y r dy 
Jo Jo Jr 

and this is finite when 2H — 2 + 2Hr — H > — 1 which gives r > jjj — \- 

The term denoted by B can be treated as follows. 
B = f 1— ^ n{n + l)(2H(2H - l)) n+1 / f duid Vj 

xe 2 e 2 l [0i „ l] (ti 2 ,---,«n+i)® n l[o ) w](wi.^.--- > «n+i)® n llh-Wi| 

z=l 

,2n /-T />T /<T >>T 



" (2ff(2fl - 1)>2 /. (T^rgt^i i i / ******** 



J1..1H „2 2£T 



xe * 2 1 e * ? l [0iUl] (n2)l [0 ,„ 2 ](wi)i?(«i,V2)" Vi-«i| 2H 2 |m 2 -v 2 | 2// 2 ■ 
Since 

n>l ^ n>0 
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we obtain 

B = (2H(2H - l)) 2 f / x 2 [ [ [ [ du x du 2 dv x dv 2 

JR. \ L + x ) JO JO JO JO 

x 2 \u 1 -v 2 \ 2H \2H-2\ \2H-2-, I m t \ 

Xe 2 |ui-t;i| \U2-V2\ l[0,«i]( u 2)l[0,^](^l)- 

We calculate first the integral du2 and dv\ and then in dx we use the change of variables 
x\u\ — v 2 1 H = y and we get 

r 1 fT f-T 2 1 ,2H 

2 / -L 2 ' ' - 



B = / jrT ^ y x J q J q du ld v 2 e 



X 

t-T 



(\ \2H-l 2H-l\ (\ \2H-1 2H-l\ 

[\ui-v 2 \ -tii ){\ui-v 2 \ -v 2 J 



= J d Ul dv 2 - ^l 2 ^" 1 " U\ H - 1 ) (\ Ul - V 2 \ 2H - 1 ~ f^" 1 ) |«i - V 2| 2 ^- 2 ^^ 

x / yV^" (|«i - w 2 | 2H + y 2 ) dy. 

JR 

which is finite if 4H — 2 + 2Hr — 3H > — 1 and this implies again r > — \- ■ 
5.2 The multidimensional fractional currents 

We will also consider the multidimensional case of the fractional Brownian motion B H in 
It is defined as a random vector B H = (B Hl , . . . , B Hd ) where B Hi are independent 
one-dimensional fractional Brownian sheets. We will assume that Hi £ (5,1) for every 
i = 1, . . . , d. In this case 

£(*) =( f 5(x- B*)dB H \. .., [ 5(x- B?)dB H A . 

\J[0,T] N J[0,T]N J 

where the above integral is a divergence (Skorohod) integral with respect to the fractional 
Brownian motion. We mention that the canonical Hilbert space Ti-H k of the fractional 
Brownian sheet B Hk is now the closure of the linear space of the indicator functions with 
respect to the inner product 



N 



<l[o,t], M°,s])n Hk = V(B? k B^) := R H *(t,s) = ]jR H ^t i ,s i 



i=i 



if t = (ti, . . . , tjv) and s = (s±, . . . , sn). In this case we have 

Proposition 4 Let B H be as above. Then for every oj the fractional Brownian current £ 
belong to the Sobolev space H~ r (R d ;M d ) for every r > max^i,. , t d — 1 + 2^) ■ 
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Proof: We write 



5(x-Bf)dBf" = [ 5 k (x-Bf)S(x k -Bf")dBf" 

J[0,T] N J[0,T] N 

= J S k (x - Bf) £ < k k (sK k (lfo3(-)) dBfK 



n k >0 



and 



n k >0 



e L [0,s] [ > 



Here again the integral dB Hk denotes the Skorohod integral with respect to the fractional 
Brownian sheet B Hk and I Hk denotes the multiple integral with respect to B Hk . Then 



X^ x 2 \s\ 2H k 



n fe >0 



A=l,l^k 



+ 



XT ^ n k 



k 

n k +l 



xl k xVgt 



sm 



n fe >0 



and 



E 



£> fc + l)! 

n fc >0 



2n fe 



KO 2 

2n fe 



cos( 2 ) * Ifo^) 



n fe >0 



sm 



'-\s\ ZH k 



P 2 1«l™fe/'. N | 



(H H , )®(»fc+i) 



|2H fe -2 



dudv 



Note that, if /, g are two regular functions of iV variables 

(/,(/)„ = (H k (2H k -l)) N f [ f(u)g(u)H\ 

J[0,T] N J[0,T] N fj[ 

:= (H k (2H k -l)) N f [ f(u)g(v)\\u-v\\ 2Hk - 2 dudv 

J[0,T] N J[Q,T] N 

if u = (ui, .., un) and v = (vi, .., vjv) and if F, G are regular functions of N(n k + 1) variables, 



fc=i 

/ / F(u 1 ,..,u nfc+1 )G(t; 1 ,..,w nfc+1 ) TT ||n i -^|| 2Hfc -W..dn nfe+1 ^ 1 ..dt; nfc+1 

J[0,T] JV ("fc + 1 ) J[0,T] JV ("fc + 1 ) 



i=l 
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and thus, by symmetrizing the above function and taking the scalar product in (%J 8 ' n,!+1 ' 
e||*(x) 

= y x * (H k (2H K - i))Mn*+i) y du dv TJ 

f^ (n fc + l)! y [0 ,r]^(» fc +i) y[o,r]^(«fc+i) ^ 



«'-« , || 2H *- 2 



x cos x i B u* cos X/ XlB vJ \ e k 2 e h 2 

\l=l,l^k J \l=l,l^k J 

x ijj* (u\ . . . , tf, . . . , « Bfc+1 )lg3] i< J r' 

+ y x k (H k {2H K - y du dv\\ 

«T>o (nfe + ^ ij^i-/[o,n JV(B * +1) ^[o,r]^(» fc +i) ^ 



I Z Jii2.ffi.-2 



sin [ x,sj sin arjBj ] e 



2 e 



Here Itz*! 2 ^ = (u\...u % N ) 2Hk and l[o jU «] = i [o,M ! 1 ]--l[o,u l Ar ] ^ u * = ( u ii--> u at)- The next step is 
to majorize |(cos(u) cos(t>)| by ^(cos 2 (u) + cos 2 (w)) and similarly for the sinus. 

e|4(x)| 2 



2n k n k +l n k +l 

< cst . y x k ( Hk (2H K - 1))^+D y du dvT] 



I I l\\2H k -2 
\U — V \\ k 



n k >0 

e~ xI ^e- X -^l^ ] (u\...,u\ 

2n k n k +l n k +l 

+csL y X k (H k {2H K - l)) N ^+V y du dvT] 

„f^0^ + 1 ) ! i^ih°n N{nk+1) J ^ N(nk+1) 1=1 

e-^^e-^^l^iu 1 , ...,u\.. .,u n " +1 )l^ i] (v 1 , . . . . . . := C k + D k 

and the two terms above can be treated as the terms A and B from the one-dimensional 
case. For example the term denotes by C k (we illustrate only this term because D k is similar 



\u l - v l \\ 2H «- 2 
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to the term B in the one-dimensional case) gives 

2n k r r n k +l 

C k = J2 X -±-[H k {2H K -l)f^ / / Y[\\u l -v l f Hk - 2 



n fc >0 "L^J * "L^J > ' z=1 

2 i l|2H t . 2 



[0,T] N< - n k+ 1 ) J [o,T] Ar ( n fc+ 1 ) f: 

e i [Qul] [u ..,u ) L [ y]\ v >•••>•••> v ) 

2n k 2 I l|2H" fc 2 I l|2H fc 

y ^(^(2^-1))^ / / e-^—e- k! ^—\\u 1 -v 1 \\ 2H *- 2 R*>'(u L ,v 1 ) n '>d 1 dv 1 

I- I- 2 1 li2H t 2 I li2fl" t . 

(H fc (2H fc - 1))" / / e^^e^^^e-^^—^ -v'f^d'dv 1 

J\0,T] N J\0,T] N 



'[0,T] N J[0,T] N 

withi^u 1 ,*; 1 ) = R^iu^vD-.R^iu^^lj). But R Hk (u,v)-u 2Hk -v 2Hk = E(B^ k B^ k )- 
E(B^ k ) 2 - B(B^ k ) 2 = -E(B^ k - B^ k ) 2 for any u,v 6 R n and this implies 

dx{l+\x\ 2 )- r C k = est. [ dx{l+\x\ 2 )- r [ f e- * fcE(B "2~ B " fc> \\ u -v\\ 2Hk - 2 dudv. 

jR d J[Q,T] N J[0,T] N 

The case N = 1 can be easily handled. Indeed, in this case 



dx(l + \x\ 2 y r C k = / / \u-v\ 2Hk - 2 I {l + \x\)- r e-^^- dx 

l[0,T] J[0,T] 



xhu-v\ 2H k 



lu _ v] 2H k -2 + 2H k r- d H k j - % ( ^ _ y] 2H k + y 2y r dy 

[0,T] J[0,T] 



»2 



which is finite if r > ^7 + f — 1- When N > 2 we will have 



/ dx(l + |x| 2 )- r C fc 

= / / ||« - v\\ 2H *- 2 (E(5*» - ^) 2 ) l+r / e-^/ 2 [E(itf* - B« k f + y 2 ]- r ^ 

./[0,Tp J[0,T] N JRd, 

and this is finite if r > ^jj- + | — 1. 

Remark 3 7n i/ie authors obtained the same regularity in the case of the pathwise 
integral with respect to the fBm. They considered only the case H k = H for every k = 1, ., d 
and N = 1. On the other hand, f^j the Hurst parameter is allowed to be lesser than one 
half, it is assumed to be in (4, 1). 



6 Regularity of stochastic currents with respect to u 
6.1 The Brownian case 

We study now, for fixed the regularity of the functional £ in the sense of Watanabe. 

We would like to see if, as for the delta Dirac function, the stochastic integral £(x) keeps 
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the same regularity in the Watanabe spaces (as a function of uj) and in the Sobolev spaces 
as a function of x. We consider B a one -dimensional Wiener process and denote 

£(x)= [ 5{x-B s )dB s 

J a 

with a > 0. Recall that by ((10]) 

ax) = J2 r n+i ((<( S )i;; s] (-)i M ( S ))] 

n>0 

(we consider the integral from a > instead of zero to avoid a singularity) where 

1 n+l 



i=l 



where as above t{ means that the variable t{ is missing. 

We recall that F is a random variable having the chaotic decomposition F 
Yln^nifn) then its Sobolev- Watanabe norm is given by 



(«)■ 



(21) 



We will get 



We use the identity 



^(n + 2r(n + l)![|(a*( a )l^ ] (.)J ll^^ 

n>0 

]T(n + 2) Q n! f L s (x)H n 



ri 4 >0 



2 o /'OO 

^n(2/)e"^ = (-l)["/ 2 l2t— / u n e~ u2 g(uyV2)du 

n!7r y 



(22) 



where g{r) = cos(r) if n is even and g(r) = sin(r) if n is odd. Since \g(r)\ < 1, we have the 
bound 



H n (y)e 



n 2 „,n + L 
<22 r :=Cri 
n!7r 2 



Then 



||£(z)||l, a < est. J>+ 2)^ /^d. 



(23) 
(24) 



and this is finite for a < -g- since by Stirling's formula n!c 2 behaves as csi.^=. 
We summarize the above discussion. 
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Proposition 5 For any x G R i/ie functional £(x) = 5(x — B s )dB s belongs to the Sobolev 
-Watanabe space D _a ' 2 /or any a > I. 

Remark 4 ^4s /or £/ie defta Dirac function, the regularity £ is i/ie same wrai/i respect to x 
and with respect to lo. 

6.2 The fractional case 

In this paragraph the driving process is a fractional Brownian motion (B^ )t£\o,T] with 
Hurst parameter H G (^,1)- We are interested to study the regularity as a functional in 
Sobolev- Watanabe spaces of 

£(z)= F 5(x-B?)dB? (25) 
J o 

when x G R is fixed. We will that now the order of regularity in the Watanabe spaces 
changes and it differs from the order of regularity of the same functional with respect to 
the variable x. 

We prove the following result. 

Proposition 6 Let (B^ )te[o,T] be a fractional Brownian motion with Hurst parameter H G 
(i, 1). For any x G R the functional (0|) is an element of the Sobolev - Watanabe space 
1D)- Q ' 2 with a > | - 

Proof: We will have in this case 

n^)iiL=E( n + 2 )> +i ) ! ii (^^s^rii^ 

n>0 H 

where we denoted by Tin the canonical Hilbert space of the fractional Brownian motion. 
As in the proof of Proposition 2 we obtain 

l®n ( \\~ i|2 



n+1 „ „ n+1 



(ff (2if - ^ f f , , , TT , ,2^-2 

7 — — tto > / / dui . . . du n+1 dvi . . . dv n+ i \\\ui-vi\ 

(n+1) 2 .j^ y [0 ,T]»+i 7[0,T]"+i £J 

xa*(ui)a^(uj)l^.](wi, • ,w n +i)l^.](vi,. . . . . . ,« n +i) 

(iZ"(2Jf - l)f +1 ^ /■ /■ , ^r 1 , l2 „_ 2 

t — — tto > / / du x ... du n+1 dv 1 . . . dv n+ i \u t - v t \ 

(fl + l) 2 ^ J[0,T]"+i 7[0,T]"+i f =1 

xal(ui)al(vi)l^ Ui] (ui, ...,u h .. . , n n+ i)l^](ui, .,u n +i) 

(tf(2# - 1))™ +1 /■ /■ J , TTl .2H-2 

H t — — tto V / / aui . . . dw n+ idt»i . . . av n +i I - vi\ 

( n + l ? i¥ ^ j=l J[0,T)n+l /[O.Tjn+l 

xa^(u i )a^(v J )l^. ] («i, ... ,u h .. .,^+1)1^.^1,.. . . . . ,v n +i) 

A(n) + B(n). 
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The first term A(n) equals, by symmetry, 

A(n) = {H{2H - [ T [ T dudva x n (u)a x n (v)\u-v\ 2H ~ 2 ( f f \v! - vf H ~ 2 du'dv' 

n + 1 Jo Jo \Jo Jo 



and using equality (|20j) we get 

A(n) = (^-D) f T f T dudv 
n + 1 Jo Jo 



x\u-v\ 2H - 2 R(u,vr U - Hn v- Hn Pu2H (x) Pv 2H(x)H n (4^) H n 



Using the identity (|22|) 



n 

n>0 

T f-T 



\u-v\ 2H - 2 



c(tf) V(n + 2) a ra! [ [ dudv\ 
n>o Jo Jo 

xR(u,v) n u- Hn v- Hn p u2 H{x) Pv 2H(x)H n H n (JL 

rT 

c(H)y^{n + 2) a n\c 2 n / dudv\u - v\ 2H ~ 2 R{u, v) n u~ Un v- Hn vr H v~ H 

JO Jo 



By the selfsimilarity of the fBm we have R(u,v) = u 2H R(1,^) and by the change of 
variables v = zu we obtain 



Y^{n + 2) a (n + l)\A{n) < 

n>0 

c(H) £(„ + 2)"nU> ( f V** du) f R{ \f (1 = dz. (26) 

n>0 V^O J Jo Z 

Using Lemma 2 (actually a slightly modification of it) in [2] 

/W^*,^ (27) 

Jo * z n2i? 

and the right hand side of (|26|) is bounded, modulo a constant, by 

y^(n + 2) a n!c^n _ 27? 

n>0 

and since n!c^ behaves as -t=, the last sum is convergent if^ — > 1, or —a > — + 
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Let us regard now the sum involving the term B(n). It will actually decide the 
regularity of the functional Following the computations contained in the proof of 

Proposition 



n>0 

rj~\ rr\ rj~i rj~\ 

= c(H) V(n + 2)"(n + l)! " (n+ i . 1 2 ) [ [ [ [ d Ul du 2 d Vl dv: 
f-f {n + 1J Jo Jo Jo Jo 

xl{o,u 1 ](u 2 )l[o,v 2 ](vi)R(u 1 ,v 2 ) n ~ 1 \u 1 - vi| 2 ^" 2 |n2 - v 2 \ 2H ~ 2 a x n {u 1 )a x n {v 2 ) 

< c{H)^{n + 2) a nn\c 2 n 

n>0 

T pT / pm \ / pv 2 

J du 1 dv 2 R{u 1 ,v 2 ) n - 1 Ui Hn V2 Hn (J \u 2 -v 2 \ 2H - 2 du 2 \ (J |m- V i| 2 ^- 2 

< c(H) V(n + 2) Q nn!c 2 f f d Ul dv 2 R{u u v 2 ) n - 1 uJ Hn v^ Hn . 

n>0 J ° J ° 

Using again Lemma 2 in [2] , we get that the integral 

f T f T du 1 dv 2 R(u 1 ,v 2 ) n - 1 u^ Hn V2 Hn < c{H)n-^H 
Jo Jo 

and since the sequence n\c\ behaves when n — > oo as we obtain that the sum X^n>o( n ~'~ 
2) a {n + \)\B(n) converges if —a — ^ + > 1 and this gives —a > | — jg. ■ 

Remark 5 Only when H = ^ we retrieve the same order of regularity of \25\) as a function 
of x and as a function ofuj. 
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